Beam debunching due to ISR-induced energy diffusion 
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One of the options for increasing longitudinal coherency of X-ray free electron lasers (XFELs) is 
their seeding with microbunched electron beam. Several schemes leading to significant amplitude of 
the beam bunching at X-ray wavelengths were recently proposed. All these schemes rely on beam 
optics having several magnetic dipoles. While the beam passes through a dipole, its energy spread 
increases due to quantum effects of synchrotron radiation. As a result, the bunching factor at small 
wavelengths reduces since electrons having different energies follow different trajectories in the bend. 
We rigorously calculate reduction in the bunching factor due to incoherent synchrotron radiation 
while the beam travels in arbitrary beamline. We apply general results to estimate reduction of 
harmonic current in common schemes proposed for XFEL seeding. 

PACS numbers: 41.60.Ap, 41.85.Ja, 42.50.Wk, 41.60.Cr, 52.59.Wd 
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I. INTRODUCTION 



Free electron laser (FEL) can be scaled to generate 
narrowband radiation in a wide frequency range. Cur- 
rently existing FELs cover bandwidths from THz to hard 
X-rays. Conventional method for decreasing the FEL 
bandwidth implies placing the undulator inside the opti- 
cal cavity. However, this approach fails for X-ray FELs 
(XFELs) due to lack of high quality reflecting mirrors at 
these frequencies. As a result, modern XFELs operate in 
a single pass self-amplified spontaneous emission (SASE) 
regime in which broadband shot noise due to discreet- 
ness of the electron beam is amplified within the FEL 
bandwidth jlj. The bandwidth of these XFELs can be 
reduced through their seeding with a narrowband signal 
rather than with a white noise. 

Growing FEL mode couples radiation with particle mo- 
tion. Therefore, XFEL can be seeded either with a nar- 
rowband radiation or with electron beam modulated 
at the X-ray wavelength Q ■ The second option looks at- 
tractive since electrons are charged particles and they 
can interact with electromagnetic fields unlike radiation 
which only weakly interacts with materials at X-ray fre- 
quencies. Several schemes predicting significant bunch- 
ing at high harmonics of available coherent light sources 
were recently proposed [B-0]- AH these schemes utilize 
magnetic bends for manipulations with modulated elec- 
tron bunch. Electrons passing through bend emit syn- 
chrotron radiation. At high energies electrons emit radi- 
ation at high frequencies and quantum effects should be 
included to find correct electron energy loss. Incoherent 
synchrotron radiation (ISR) has nearly 100% frequency 
spread which results in the electron energy diffusion when 
quantum effects are accounted since electron energy loss 
can be described within random walk model. Electron 
energy diffusion translates into diffusion along longitu- 
dinal position since bends are dispersive elements (elec- 
trons with different energies travel along different trajec- 
tories). As a result, prebunched electron beam smears 
out when it passes through bends and the bunching fac- 
tor reduces. 



The effect of the beam debunching becomes stronger 
for FELs seeded at smaller wavelengths. First of all, en- 
ergetic electron bunch is required to generate short wave- 
length FEL radiation which increases ISR-induced energy 
diffusion in the beamline bends. At the same time, the 
associated longitudinal spread should be smaller to cause 
significant debunching of shorter wavelength modulation. 
These two effects combined result in stronger smearing of 
harmonic current of short wavelength XFELs. Therefore, 
there is a technological limit on the shortest wavelength 
bunching which can be created in various schemes for 
generating high harmonic content. This effect was esti- 
mated for various XFEL seeding schemes [1, H, but 
these estimates remain qualitative or numerical and do 
not allow for detailed trade-off studies. It is the purpose 
of this paper to study this effect rigorously and find ac- 
curate quantitative estimate on degradation of harmonic 
current in various bcamlines. 



II. QUALITATIVE ESTIMATE FOR SMEARING 
OF HARMONIC CURRENT 



First, we estimate the parameters region in which 
smearing of harmonic current due to ISR-induced en- 
ergy spread is significant. We consider any FEL seeding 
scheme for creating longitudinal bunching. Consider the 
last dipole of the beamline optics which recovers imposed 
modulation as longitudinal bunching. 

Quantum effects in ISR at high energies result in en- 
ergy diffusion as described in Ref. floj 
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where D is the energy diffusion coefficient, p = 
P"/mc/ (eB) is the electron gyroradius in the magnetic 
field B; h is the reduced Planck constant, and s is the 
path length. Two identical electrons entering the dipole 
exit at different longitudinal positions since their ener- 
gies become different inside the dipole due to ISR and 
the dipole energy dispersion transforms this energy dif- 
ference into difference in longitudinal positions. As a 
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result, bunched beam smears out. We approximate this 
effect assuming energies of two electrons to be constant 
along the bend but to be different by th e overall induced 
energy spread due to ISR, AE = \/2Ds. The path length 
of an electron inside the magnet of length L is equal to 



s = p arcsm — . 

P 



(2) 



Then the path difference of electrons having different en- 
ergies is 



As 



ds dp A ^ L 2 AE 



(3) 



This path difference results in smearing of the imposed 
modulation. Smearing is significant if the path difference 
is on the order of half wavelength of modulation. Then 
one obtains the qualitative estimate for parameters re- 
gion in which ISR-induced debunching is not significant. 



55 



he 2 a 7 7 5 



M^m?<? A 2 90 ' 6 A 2[^ [lQGev) 



al[deg] ( E V 



«1. (4) 



Presented estimate shows that beam bunching for hard 

X-ray FELs (A ~ 1 A, E ~ lOGcV) requires the use of 
optical elements with weak bends having sub-degree bend 
angles. However, this estimate does not predict how fast 
the bunching degrades along the bend. In particular, it 
is not clear how fast the bunching drops to zero and how 
far one can go beyond this limit without dramatic loss in 
bunching amplitude. 



III. SMEARING OF HARMONIC CURRENT IN 
ARBITRARY BEAMLINE 

A. Vlasov equation 

Any electron beam can be described as an ensemble of 
electrons occupying some phase space volume. This en- 
semble can be described with the phase space distribution 
function /(C), where £ is the 6D phase space coordinate 
of each electron 



C = (x,Px,y,P y ,At,-AE), 



(5) 



where x and y are the transverse electron coordinates in 
respect to the reference trajectory, p x and p y are the cor- 
responding momenta, At is the deviation of arrival time 
to position s along the beamline, AE is the deviation of 
particle energy from the average bunch energy. 

The bunch interacts with electro-magnetic fields while 
it travels along the beamline. These forces satisfy 
Maxwell equations which indicates that beam dynamics 
is Hamiltonian and it is fully described with Hamilto- 
nian -ff (£, s). The evolution of the distribution function 
satisfies Vlasov equation which can be considered as a 
continuity equation in the phase space 



ds 



dsf(C, s) + {/, H} = d s f + (V/) T J{VH) = 0, (6) 



where {f,H} = (V/) T J(ViJ) is the Poison bracket, J 
is the unit block-diagonal antisymmetric symplectic ma- 
trix, V is the 6D gradient in the phase space, and super- 
script T stands for transposition. 

Typically the electron bunch can be considered well lo- 
calized and quasi-monoenergetic. Under this assumption, 
forces can be well approximated to be linear function in 
respect to the phase space coordinates. The correspond- 
ing Hamiltonian describing beam dynamics is quadratic 
then, H(C,s) = l/2C T H(s)C, U = U T . The trajectory 
of each electron defines the map in the phase space which 
can be described with the transform matrix R(s, sq) for 
linear beamlines. Then formal solution for Vlasov equa- 
tion © can be presented in the following form 



f(C,s) = f(R- 1 (s,s )C,s 
dR(s,s ) 



ds 



jy.R(s, s ), R(s ,s ) = I. 



(7) 

(8) 



Alternatively, the dynamics of modulated beams can 
be conveniently described in the spectral domain as il- 
lustrated in Ref. [llT | . This representation is particularly 
useful for description of modulated beams since they are 
well localized in the spectral domain in case of quasi- 
monochromatic modulations. The beam can be fully de- 
scribed with its 6D spectral distribution which evolution 
is Hamiltonian in case of linear optics 



; (k, S ) = J f(C, S )e lkT<: d 6 C, 



(9) 



^ = ^ + {/ k) if k } = 0, ff k = -VjWJk.(10) 

ds ds Z 

The solution of the spectral Vlasov equation was found 
in Ref. [ll| 



/ k (k, so) = / k (i? T (s,s )k,s ). 



(11) 



This solution indicates that the spectral distribution 
function remains constant along characteristics in the 
spectral domain 



k(s) = ir T (s,s )k(s ), 



(12) 



which reduces evolution of each spectral component to 
linear transform of its modulation wavevector. 



B. Boltzmann equation in phase space domain 

Synchrotron radiation of highly rclativistic electrons is 
confined within a small angle in respect to the electron 
instantaneous velocity, A9 ~ 1/7 <c 1. Therefore, emis- 
sion of a photon with angular frequency lo mainly results 
in the reduction of the electron energy by AE = hoj. 
The change of the electron transverse momenta is on the 
order of I/7 smaller then the change of the longitudinal 
momentum and can be ignored in the first order. Assum- 
ing that the emission process is instantaneous in time, 
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one can use Vlasov equation © to describe electron dy- 
namics between photon emission events. Electron energy 
loss due to photon emission results in the collision oper- 
ator which changes phase space density. Large number 
of electrons within the bunch and ergodicity property al- 
lows one to reduce discrete emission events to continuous 
change of the ensemble distribution. As a result, evo- 
lution of the electron bunch can be described with the 
following Boltzmann equation 



chrotron radiation power spectrum. 



'1 

ds 



-f = d s f + {f,H} = C^{f} 



(13) 
(14) 



dui 



where C ISR [/] is the collision operator for ISR energy 
loss which describes detailed balance between electron 
states. H is the beamline Hamiltonian if synchrotron 
radiation is neglected, (dN / dui)dujds is the probability 
for a photon emission within bandwidth duj while electron 
travels distance ds. The probability of photon emission is 
related to the power spectrum dPj duj of a single electron 
synchrotron radiation [l2j 



1 v^e 2 ui 
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dN _ 1 dP 
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K 5/3 (x)dx, (15) 



(16) 



Note that collision operator C [/] describing ISR- 
induced energy losses is linear in respect to the electron 
distribution function. This property reflects the fact that 
ISR is a single particle effect. Superposition of radiation 
fields from many particles and influence of the resulting 
Coherent Synchrotron Radiation (CSR) force on parti- 
cle dynamics is not included in this model. Presented 
model also assumes that photons do not interact with 
other electrons once they are emitted. 

Boltzmann equation described by Eq. (fTB"]) can be used 
to determine the main parameters of ISR. For example, 
one can easily find the overall energy loss and the energy 
diffusion 

^P) = [EC^{f}d e C = -P, (17) 

ds y ISR j 

(^P) = f(E-(E)?C™{f}d e C = 2D,(m) 

\ ) ISR J 

where (E) = j Efd e £ is the mean electron energy and 
(Ai? 2 ) = j(E — (E)) 2 fd 6 C is the rms energy spread. 
These coefficients can be expressed through the syn- 
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D = 



dP 

du 



-dui 
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Note that the overall change of beam energy and en- 
ergy spread are described with Eqns. (fl~7f and (fl~8|) if the 
beam is not accelerated and the beamline does not con- 
tain elements having nonzero Rqi, i = 1..6 beam matrix 
elements, e.g high order mode RF cavities introducing 
energy slews. 



Boltzmann equation and its solution in spectral 
domain 



Boltzmann equation described with Eq. (fTB")) — ([14")) 
is a complicated integro-differential equation. On the 
other hand, this equation take much simpler form in 
the spectral domain. The appropriate equation can be 
found by taking Fourier transform of Boltzmann equa- 
tion (fl3l) . Its left-hand side is the same as Vlasov equa- 
tion © and, therefore, its Fourier transform reduces to 
spectral Vlasov equation (|10p . The Fourier transform of 
the ISR collision operator results in the corresponding 
collision operator in the spectral domain. One can find 
it to be equal to 



ISR 



f(C)e 



l)^L d ojd 6 C. (21) 
duj 



These integrals can be evaluated under assumption 
that the electron bunch is quasi-monoenergetic. In this 
case the ISR spectrum can be assumed to be the same for 
all the particles and the integral over photon frequency 
can be evaluated independently from the integral over the 
phase space which yields to the spectral distribution func- 
tion. This assumption is equivalent to neglecting with the 
radiation cooling effect which is a valid approximation for 
high brightness linear machines. 

The integral over photon frequencies can be evaluated 
using spectral distribution of the ISR power (fT5j) and the 
spectral Boltzmann equation takes the following form 



ds 

ISR 



P cS (x) 
D eS (x) 



ds 



+ {h,H k } = Cl SR , (22) 

4^ cS )/k, (23) 
54 5 WTT>? - 3 sinh(|asinh(xr)) 

P, (24) 

D, (25) 

h E huj c . (26) 
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One can note that Boltzmann equation (|22[) — (|23l) in 
the spectral domain has much simpler form compared 
to Boltzmann equation in the phase space domain (fT3|) 
- f| 14[) . ISR collision operator in the spectral domain 
C k SR algebraically depends on the spectral distribution 
function which indicates that different spectral compo- 
nents evolve independently from each other. The col- 
lision operator depends on two parameters, namely the 
effective energy loss P eS and diffusion D eS coefficients, as 
described by Eqs. (|2~4"]) — (f2l)]) . The effective energy loss 
term results in the change of each spectral component 
phase which corresponds to the loss of the mean bunch 
energy while absolute amplitude of each component is not 
affected. On the contrary, the effective energy diffusion 
term results in the reduction of each spectral component 
amplitude which manifests as smearing of harmonic cur- 
rent in the phase space domain. 




FIG. 1: (Color online) Dependence of effective energy loss 



The effective energy loss and diffusion parameters de- 
pend on a single parameter x = kEhuj c . The nonzero 
energy modulation wavenumber Ue indicates presence 
of energy bands in the phase space distribution sepa- 
rated by AEbands — 2n/kE- Electrons passing through 
the bend emit photons which reduces electron energy 
by Ai?igR = fa in a single emission event. If all elec- 
trons would lose the same energy in a single photon emis- 
sion event (AE'isr) = hw c then the distribution function 
would not change if (Ai?igri) = AEt, an d s since all the 
electrons would move into the same position within the 
next band and all the bands are equally populated. Then 
energy diffusion would vanish at kE^ c = 27r in this hy- 
pothetical scenario. However, photons emitted due to 
ISR have nearly 100% frequency distribution, so in real- 
ity electrons do not move into the same position within 
the next energy band but also acquire some random mis- 
placement which manifests as energy diffusion. At the 
same time, the diffusion coefficient at >c — lix is much 
smaller compared to the case of x < 1 (as illustrated 



in Fig. Q} since a large fraction of electrons do not sig- 
nificantly change their relative position within the corre- 
sponding band. 

Boltzmann equation in the spectral domain is a linear 
hyperbolic equation which can be solved using method 
of characteristics. The characteristics describing trajec- 
tories in the spectral domain are the same as for Vlasov 
equation (|12[) . This property indicates that transform 
of modulation wavevector is independent from whether 
ISR is accounted or not. Taking ISR into account (both 
overall energy loss and energy diffusion) results in the 
change of the harmonic current amplitude and phase 
along the beamline. At the same time, the optics required 
to transform initial modulation into required one and the 
bandwidth of that modulation are not affected by ISR. 
This property is highly useful in studies of XFEL seeding 
schemes since the beamline optics can be designed with- 
out accounting for ISR which can be included at the last 
stage while estimating reduction of the output bunching 
factor. 

The absolute amplitude of the spectral distribution 
changes due to effective energy diffusion effect but it is 
not affected by the overall energy loss. One can find so- 
lution of linear Eqs. (|22j) — ([23]) in the following form 
which can be used to find smearing of harmonic current 
due to ISR in any arbitrary linear beamline 



/ k (k(s/),s/) 



/k(k(s ),so) 



cxp 



k E (s)D eS ( S )d S I , (27) 



R T (s,s )k(s )- 



(28) 



There are several properties which follow from Eq. (|27|) . 
First of all, the harmonic current amplitude always re- 
duces along the beamline. This property reflects the dif- 
fusive nature of ISR which smears out small scale vari- 
ations in the distribution function. Another property 
states that the rate of modulation smearing is propor- 
tional to k E . This property reflects the fact that ISR- 
induced diffusion is the energy diffusion. Uniform energy 
distribution is not affected by this effect. Also note that 
the integral along the beamline in Eq. (|27j) can be sub- 
stituted as a sum of integrals along each element. As a 
result, the overall attenuation of the harmonic current is 
equal to the product of attenuations in each dipole. 



D. Fokker-Planck approximation 

Boltzmann equation in the phase space domain (|13|) - 
(|14l) can be simplified under assumption that the charac- 
teristic energy scale of the distribution function / /OeI 
is much larger than the characteristic energy of emitted 
photons Hlu c . Then the distribution function in the col- 
lision operator (fT4"]l can be expanded and the following 
Fokker-Planck equation recovered 



'1 

ds 



d E (Pf) + d 2 EE (Df), 



(29) 
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where energy loss P and diffusion D coefficients are de- 
scribed with Eqns. (JT7J) and (f!5]l. 

Fokker-Planck equation can be solved in the spectral 
domain in same way as it was done in Sec. lIII Cl for Boltz- 
mann equation. However, this analysis is not required 
since Fokker-Planck equation is an approximation of a 
more general Boltzmann equation and its solution can 
be recovered from a general solution. Boltzmann equa- 
tion (fT3|) yields to Fokker-Planck approximation when 
the typical energy of emitted photons is much smaller 
than the characteristic energy scale of the electron dis- 
tribution function. This assumption corresponds to the 
limit of fc^Swc <C 1 as discussed in Sec. MI Cl Therefore, 
solution of Fokker-Planck equation (|29p can be found as 
a small photon energy limit of the general solution for 
Boltzmann equation (|27p. when D cS = D. 

Smearing of harmonic current depends on the energy 
modulation wavenumber along the beamline which can 
be found using transform (|28[) . Consider, for example, 
some beamline for XFEL seeding scheme. As discussed 
in Ref. [Til ] such a beamline is designed to recover im- 
posed modulation as longitudinal bunching at a given 
wavenumber, k(s/) = k 2 fc = (0, 0, 0, 0, k, 0) T . Then 
the modulation wavevector at any given position along 
the beamline can be expressed through the modulation 
wavevector at the final position and the linear transform 
matrix as k(s) = R T (sf, s)k(s/). Then one can find 
that the overall attenuation of harmonic current due to 
ISR-induced energy spread in the Fokker-Planck approx- 
imation is equal to 



/k(k(s/),s/) 



/k(k(s ),s ) 



exp 



-k 2 [ DR 2 56 (s f , 



s)ds 



(30) 

This result agrees with qualitative estimate presented in 
Ref. 0. The solution shows that the degradation of 
harmonic current strongly increases in beamlincs having 
large energy dispersion. Therefore, conventional XFEL 
seeding schemes utilizing chicanes might be affected 
by ISR since chicanes are specifically designed to have 
large energy dispersion R$q. 



tion wavevector in the spectral domain. The trajectory 
of each spectral component does not depend whether ISR 
energy diffusion is neglected or taken into account which 
follows from Eq. (j2"8")l since it is independent on the diffu- 
sion coefficient. Taking into account ISR-induced energy 
diffusion manifests as the next order effect and results 
in reduction of the harmonic current along the beamline. 
Then the exponential factors in Eqs. (|2"T|) and (|3"0"|) can 
be interpreted as attenuation of harmonic current due to 
ISR-induced energy diffusion. 

In our analysis we limit ourselves to modulations which 
are recovered as longitudinal bunching at required wave- 
length at the end of a chosen beamline element. Then 
the beam modulation at the beginning of the beamline 
element is well defined since its transform is described 
with Eq. (|28[) . From this prospective, the precise mech- 
anism of imposing initial modulation is not important. 
Therefore, the attenuation factor for each element will 
depend on the output bunching wavelength and beam- 
line parameters. 

We will present estimates for the attenuation factor 
based on the solution of Fokker-Planck equation de- 
scribed with Eq. (|30[) since it is much simpler than 
solution of Boltzmann equation (|27|) and the integrals 
can be calculated analytically. This approximation is 
valid when condition fcs/iw c < 1 is held. The energy 
wavenumber of modulation can be expressed through 
the wavelength of final longitudinal bunching and the 
beamline dispersion. In this section we will use conven- 
tional units used in Beam Physics, particularly we will 
use relative energy spread instead of full energy devia- 
tion, AE = -fmc 2 (A-f/-f). Then condition for validity of 
Fokker-Planck equation expressed in common units reads 
as 

(R 5 6(s f , s)) max « ^3 = 2.39 p L J, - urn, (31) 
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B[T] E 



where A is the output bunching wavelength and A e 
2irh/(mc) ps 2.43 • 10~ 12 m is Compton wavelength. 



A. Single bend 



IV. SMEARING OF BUNCHING IN VARIOUS 
XFEL SEEDING SCHEMES 

In this section we calculate attenuation of harmonic 
current in various beamlincs designed for XFELs seeding 
with microbunched electron beams. In these schemes the 
beam is typically modulated in the wiggler by interacting 
with external laser at resonant wavelength and modula- 
tion at high harmonic is recovered in the following linear 
beamline as longitudinal bunching. This dynamics can 
be conveniently described in the spectral domain as was 
illustrated in Ref. This formalism describes change 
of the spectral distribution function along characteristics 
which can be considered as trajectories of the modula- 



First, we estimate attenuation of harmonic current in a 
single bend which can be the last bend of a more compli- 
cated linear beamline. The transform matrix of a bend is 
described with Eq. (|A10[) . Then attenuation of harmonic 
current can be found from Eq. (|30[) and for a single bend 
it is equal to 



ibend 



exp 



exp 



exp 



-k 2 ° J R 2 5e (0)d(p 
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fine 



m 



f E \ 5 a 7 [deg] 
VlOGcV; A 2 [A ] 



(32) 
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where a is the dipole bend angle, a fine = e 2 /(ftc) ~ 
1/137 is the fine structure constant, and r e = e 2 /(mc 2 ) w 
2.818 • 10 _15 to is the classical electron radius. 

Smearing of harmonic current can be neglected if the 
exponent argument is smaller than unity. That results in 
a maximum angle of the bend which can be used in the 
beamline. Exceeding this critical value results in a strong 
smearing of harmonic current due to very strong scaling 
of attenuation versus bend angle. The parameters region 
for importance of ISR debunching agrees with qualitative 
estimate (U). Qualitative estimate shows the same scaling 
versus main parameters as in a rigorous analysis but the 
numerical factor is about a factor of 6 larger since the 
ISR-induced energy spread was assumed to be acquired 
in the beginning of the bend rather than being uniformly 
distributed along its length. 



B. Chicane 

Most of the beam based XFEL seeding schemes uti- 
lize chicanes 0-0] • To estimate the effect of ISR in these 
schemes we consider that the chicane is designed to trans- 
form imposed space-energy modulation as final longitu- 
dinal bunching. The detailed analysis of this setup is 
presented in Appendix Eland it yields to the attenuation 
factor of 



scheme due to ISR can be estimated as 



A 



ch 



exp 




6700 



a 7 [deg] / E 



A 2 [A] VlOGeV 
0.072] +0.022 



(33) 



One can note that attenuation of harmonic current is 
stronger for larger chicanes. This result follows directly 
from the general expression ([30]) for the attenuation fac- 
tor in a general beamline. Typically, the chicane strength 
scales as R^q ~ pa 3 . Then this term becomes dominant 
compared to other numerical factors and the following 
approximate expression for attenuation can be used most 
of the time. 



.4, 



ch 



exp 



a[dcg]i? 2 6 [ M m]B 2 [T] / E 



4.7A 2 [A] 



VlOGeV 



(34) 

This expression allows one to estimate the effect of ISR in 
proposed schemes for XFEL seeding with microbunched 
beam. 

HGHG scheme Q uses a single modulator and a 
single chicane. The chicane strength required to 
achieve maximum bunching is approximately equal to 
R^i^mod/i) = MniV 2 7r, where fx n i — n is the first 
maximum of the n-th order Besscl function, J' n ([i n i) = 
II ll . Then attenuation of harmonic current in HGHG 



Ahghg ~ exp <^ -5.4 ■ 10 >nl 



( E 



\10GeV 



A 7rii 



a[deg]B 2 [T] 



Mmod) (A 7rms / 7 [0.01%]) 2 



(35) 



Note that degradation of harmonics current does not de- 
pend on the wavelength of the output bunching directly. 
It depends only on the harmonic number used in HGHG 
scheme since ~ n at n > 1. 

Similar estimates can be made for harmonic current 
smearing in EEHG @ and CHG @ schemes. These 
schemes utilize two chicanes of different strengths, so the 
attenuation should be estimated for both chicanes. The 
first chicane in both schemes is used to create energy 
bands in the phase space from the modulation which can 
be considered as mostly longitudinal bunching, kE(sf) S> 
fce(so)- Therefore, the estimate (|3"3"]l can be applied to 
the first chicane even though it was derived for a chicane 
which converts spatio-energetic modulation into purely 
longitudinal bunching. Straightforward algebra shows 
that the ratio of chicane strength and modulation wave- 
length is almost the same for both chicanes in EEHG and 
CHG schemes, (i? 5 6/A) c hicanci ~ (-fiWA) c hicanc2- There- 
fore, both chicanes in those schemes result in the same 
attenuation of harmonic current if they are designed us- 
ing the same dipoles and the difference in their strengths 
is caused by different drift lengths within the doglegs. 
Then attenuation of harmonic current in EEHG scheme 
can be estimated as 



Aeehg « exp | -10.8- 10 3 /4i 
A 7rm 



E 



\10GeV 
a[deg}B 2 [T] 



A 7mod2 ; (A 7rms / 7 [0.01%]) 2 



(36) 



where A 7mo( j2 is the energy modulation imposed in the 
second modulator. In this estimate we neglected with en- 
ergy diffusion due to ISR in the second undulator which 
modulates the beam. This effect can be easily accounted 
under assumption that the undulator energy dispersion 
is much smaller compared to chicanes. Then the en- 
ergy modulation wavenumber fc^ = /ini/(A 7mo d2/ 7 ) is 
constant along undulator and full attenuation factor for 
EEHG scheme should be reduced by exp(—k E D un dL), 
where L is the undulator length and D un d is the undula- 
tor energy diffusion coefficient found in Ref. [l~3| (the rms 
energy change by the undulator is (A 7 / 7 ) 2 = 2D un dL). 
Similar estimate for CHG seeding scheme yields to 



Achg ~ exp 



er 2 [^m] / A 7 



-4.3- 10 7 (M + 1) 2 ( 



— V 

lOGeV J 



«[dcg]B 2 [T] 



A 2 [A ] V A 7md J (A 7rms / 7 [0.01%]) 2 j ' 



(37) 
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where M is the compression factor in CHG scheme and 
Aji n d is the energy slew imposed on the beam inside 
cavity (additional energy at the location of rms bunch 
length), and a z is the bunch length after compression. 
The estimate for smearing of harmonic current in CHG 
scheme indicates that this effect is very strong compared 
to other seeding schemes. This result comes from the fact 
that chicanes used in CHG are much stronger than what 
is required for HGHG and EEHG schemes. Chicanes in 
CHG should change relative particle position on the order 
of the pulse length to provide significant compression un- 
like harmonic generation schemes which require change of 
electron relative position on the order of the modulation 
wavelength. Therefore, CHG scheme may be not feasible 
for seeding XFELs due to strong ISR-induced smearing 
of harmonic current. 



C. Emittance Exchanger 

Smearing of harmonic current is large in beamlines 
having large energy dispersion i? 56 as discussed in 
Sec. IIII Dl ISR-induced attenuation can be reduced by 
specifically designing beamline without energy disper- 
sion. This approach implies that conventional schemes 
for XFEL seeding such as HGHG, EEHG, and CHG can- 
not be implemented since they rely on the presence of 
dispersive elements to recover imposed spatio-energetic 
modulation as longitudinal bunching. Therefore, one 
should consider imposing modulation different from mod- 
ulation in E — z phase plane. For example, one can con- 
sider using Emittance EXchanger (EEX) optics J14J to 
transform transverse modulation into bunching J7|, 1 1 51 ] - 
Properly designed EEX swaps longitudinal and trans- 
verse phase spaces of the beam, and therefore, the energy 
dispersion of the beamline i?56(s/,sn) is zero. However, 
the energy dispersion from the middle of the beamline 
to the end, R^{sf,s), is not zero, and attenuation of 
harmonic current does not vanish in this setup as fol- 
lows from Eq. (|3"0|) . Additional minimization of the ISR- 
induced debunching can be achieved by using doglegs 
having zero dispersion so that energy modulation is kept 
small in the middle of the EEX optics. Zero energy 
dispersion of the dogleg can be achieved by inserting a 
quadrupolcs triplet in the dogleg drifts so that the entire 
drift looks like an effective neg ative drift space in terms 
of linear transform matrix 



1G] 



We design the following setup to recover longitudinal 
bunching using EEX optics. First, the beam is transver- 
sally modulated when it passes through the mask. Then 
the imposed modulation is recovered as longitudinal 
bunching in the following EEX optics. The dispersion of 
each dogleg is chosen to be zero to minimize ISR-induced 
smearing of harmonic current. Such a design requires 
that the beam is modulated along x' phase space coor- 
dinate in front of the EEX optics. This can be achieved 
by inserting a focusing lens just in front of the EEX and 
placing transverse mask in the focus plane of this lens 



as discussed in Ref. [16j. The attenuation of harmonic 
current in this scheme can be calculated using the same 
algorithm as discussed in Appendix [S] and one can find 



^4eex = exp 



A 2 [A] 



(38) 



We numerically verify the quantitative estimate (|38j) 
by simulating EEX optics with a simple particle push- 
ing code. All the beamline elements were considered to 
be linear and the ISR was included as a random change 
of the electron energy at each time step within bends. 
The results arc presented in Fig. [5] for the regime of 
ItEtiuJc C 1 in which Fokkcr-Planck approximation can 
be used. The following parameters of the EEX optics 
were used: beam with mean energy of E — \2GeV 
and rms energy spread of A7/7 = 1CP 4 , normalized 
transverse emittance e n = 0.14 mm-mrad, beam sizes 
a x = o v = a z = 100/im in all dimensions, EEX trans- 
forms initial transverse modulation with 23nm wave- 

o 

length into 0.3A longitudinal bunching, EEX dipolcs with 
B = IT magnetic field bend the beam to a = 0.3 degree 
angle. Fig.[2]shows perfect agreement between numerical 
simulations and quantitative analysis. 




600 800 
distance [mm] 



1400 



FIG. 2: (Color online) Attenuation of harmonic current in 
EEX optics calculated numerically (solid blue line) and ana- 
lytically (dotted red line). 

Note that attenuation of harmonic current in EEX op- 
tics described with Eq. (|38p is similar to attenuation in a 
chicane having zero energy dispersion R^e = described 
with Eq. (|33| . This fact can be explained by similarity of 
these two optics which require the same doglegs having 
zero dispersion. Moreover, the scaling of attenuation is 
the same for both schemes, \nA cx — a 7 j 5 /X 2 , but the 
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corresponding numerical factors arc different. The nu- 
merical factor for attenuation of harmonic current in the 
chicane is on the same order as for EEX if = and 
it is an order of magnitude larger compared to EEX if 
i?56 ~ pa 3 . Also note that the same scaling holds for a 
single bend described by Eq. (|32|) . 



V. DISCUSSION 

We developed quantitative approach which allows one 
to calculate smearing of harmonic current due to ISR- 
induced energy spread in an arbitrary beamline. The 
approach is based on a beam representation in the spec- 
tral domain where modulation is represented as a well- 
localized distribution which evolves along the trajecto- 
ries. The energy diffusion manifests in the spectral do- 
main as a damping operator which allows one easily cal- 
culate attenuation of harmonic current amplitude along 
the beamline. The analysis can be further simplified us- 
ing Fokker-Planck approximation which is valid when the 
typical energy of emitted photons is much smaller than 
the energy bands of modulation. 

We applied developed formalism to estimate smear- 
ing of harmonic current in various schemes proposed for 
XFEL seeding in which imposed modulation is recovered 
as longitudinal bunching in the following beamline. We 
demonstrated that attenuation of harmonic current am- 
plitude increases in beamlines having large energy dis- 
persion i?56- At the same time, considered schemes re- 
sulted in the same scaling for attenuation of harmonic 
current, InA cx -a 7 ~f 5 /\ 2 as follows from Eqs. ([32]). (f33| . 
and (f3"51) . The only difference between beamlines came 
as different numerical factor in front of the scaling. This 
scaling indicates very rapid increase of ISR-induced de- 
bunching effect when cither beam energy or the dipolc 
bend angle increases. As a result, the proposed XFEL 
seeding schemes should utilize elements with small bend 
angles. The maximum bend angle can be estimated from 
the condition that the ISR-induced diffusion does not re- 
duce the amplitude of harmonic current by more than 
a factor of 2. Rapid scaling of the attenuation versus 
bend angle indicates that different beamlines having dif- 
ferent numerical factors for attenuation scaling result in 
similar critical angles for bends. From this prospective, 
significant complication of the beamline optics to reduce 
ISR-induced debunching does not seem practical. 
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Appendix A: Calculation of harmonic current 
reduction in chicane 



In this Appendix we present algorithm for calculating 
attenuation of harmonic current in chicane which trans- 
forms initial beam modulation into longitudinal bunch- 
ing. However, the same algorithm can be used to calcu- 
late smearing of an arbitrary short-scale modulation in 
an arbitrary beamline. 

Degradation of harmonic current can be described with 
Fokker-Planck equation ((29)) which describes beam en- 
ergy diffusion coupled with the 6D phase space trans- 
port. In Sees. MI CI and MIDI we assumed phase space 
variables to be canonical conjugate. However, the vari- 
ables do not have to be canonical in order to correctly 
describe the beam transport since in our analysis of ISR 
we did not use condition for variables conjugancc. There- 
fore, one can choose conventional set of variables, i.e. 
C{s) — (x, x' , y, y', cAt, A7/7), where x' and y' are the 
particle angles in respect to reference trajectory. 

We describe the chicane as two doglegs separated by a 
drift space as illustrated in Fig. [3] To simplify our anal- 
ysis, we consider hard edge bends. Then particle motion 
in y-plane is decoupled from motion in other phase space 
planes and it is represented by a simple drift. Therefore, 
energy diffusion affects motion only in x and z phase 
space planes and the system dynamics can be adequately 
described in 4D phase space. 



— d ^S,^ d - — S 2 —» d ^S,^ d +- 




FIG. 3: (Color online) Schematics of a chicane. 

The transform matrix of the chicane shown in Fig. [3] 
is equal to 
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Here the dogleg parameters can be found in the ultra 
relativistic limit as 



L = 



2d cos a + Si 



(A3) 



5*1 + 2d cos a - (2d + Si) cos 2 a , . A . 

V = „ ; 9 : , (A4) 



sin a cos' a 
2d sin a cos a — 2da cos 2 a + Si sin 3 a 
sin a cos 2 a 



(A5) 



where a is the bend angle of the dipole with parallel pole 
faces, d is the dipole length, Si is the distance between 
bends in each dogleg, and S2 is the distance between two 
doglegs. 

Finding attenuation of harmonic current requires 
knowledge of the modulation wavevector inside each 
bend. First we find the modulation wavevector at the 
edge of each dipole using Eq. (|28j) which describes its 
transform along the beamline. We consider the out- 
put modulation of the beam to be longitudinal bunch- 
ing, k(4) = 2tt/A x [0,0,1,0] T . Then the modulation 
wavevector at positions 1, 2, and 3 along the beamline 
can be found as 



k(l) 
k(2) 



^ie g k(l) = 2^0,-77,1,^ 



i? c J hicano k(4) = 2tt/A[0, 0, 1, 2£] J , (A6) 

(A7) 
(A8) 
(A9) 



k(3) = i^ t k(2) = 27^0,-77,1,^ 
k(4) = 2tt/A[0,0, 1,0] T . 



The wavenumber of modulation inside each bend can 
be found using transform matrix of the bend from its 
edge to some intermediate plane 



Rg — 



cos 9 p sin I 

sm9/p cos 9 

— sin 9 p(cos 9 - 





— cos( 
sin0 

1) 1 p(sm9-( 
1 



where the intermediate angle inside the bend is defined 
as 9 to distinguish it from the overall bend angle a of the 
dipolcs. Then the modulation wavevector in each bend 



k«(0) = iV r k(l), 
k^(a-9) = (R mp R e R mp ) T k(2), 
k^(9) = (R mp R e R mp )- T k(3), 
kW(a-0) = i#k(4). 



(All) 
(A12) 
(A13) 
(A14) 



Using Eq. (|27| in the limit of ^Suc -C 1 one can 
find attenuation of the harmonic current in each dipole. 
Note that we use of A7/7 instead of AE as an inde- 
pendent energy variable and the bend angle 9 instead 
of s as the path length variable in Fokker-Planck equa- 
tion. This choice of variables changes the diffusion coef- 
ficient to Dg = p I '(7771c 2 ) 2 D and the energy modulation 
wavenumber to £^7/7 = jmc 2 ks- The overall reduc- 
tion in the modulation amplitude is equal to the product 
of attenuation factors in each element. Calculating ap- 
pearing integrals and considering the limit of small bend 
angles, a<l, one obtains the following expression for 
attenuation of harmonic current in the chicane 



chii 



cx p{~E/ (k { £ ih {9)) 2 Dgd9 



r [deg] / E V 



0.072 + 0.022 



pa° 



(A10) where -R56 = 2£ is the chicane strength. 



(A15) 
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